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Abstract. A notion of rational Baker- Akhiezer (BA) function related to a configuration 
of hyperplanes in C" is introduced. It is proved that BA function exists only for very spe- 
cial configurations (locus configurations), which satisfy certain overdetermined algebraic 
O ' system. The BA functions satisfy some algebraically integrable Schrodinger equations. 
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so any locus configuration determines such an equation. Some results towards the clas- 
sification of all locus configurations are presented. This theory is applied to the famous 



P_|, Hadamard's problem of description of all hyperbolic equations satisfying Huygens' Princi- 
ple. We show that in a certain class all such equations are related to locus configurations 
and the corresponding fundamental solutions can be constructed explicitly from the BA 
functions. 



^ ■ Introduction 



The notion of Baker-Akhiezer function (BA function) has been introduced by Krichever 
in the theory of finite-gap or algebro-geometric solutions of the nonlinear PDE's, inte- 
grable by the inverse scattering method . The BA function is a far-going generalisation 
of the classical function 

^ ^ a{x-z) ^^(^^^ 
well-known as a solution to the classical Lame equation: 

L^ = A^, L = -— + 2p{x), \ = -p{z). 
Here a,( and p are classical Weierstrass elliptic functions (see e.g. H). 
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In the degenerate case one has the corresponding trigonometric and rational versions: 
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^.., = (l-^cotx)e^ L = -— + -^, 
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kx ' dx'^ 



Certain multidimensional versions of these functions in the rational and trigonometric 
cases have been introduced by Chalykh and Veselov in [Q in the theory of quantum 
Calogero-Moser problem. In this paper we will restrict ourselves by the rational case only. 
The construction of (see also 0) relates such a BA function ip to a. configuration 21 of 
the hyperplanes in a complex Euclidean space C" given by the equations (a, x) = , 
taken with some multiplicities rria G Z_|_ . Here a & A, A is a finite set of noncoUinear 
vectors. The function il){k^x) , k,x E C" is determined by certain analytic properties in 
k (see Section 1) and exists only for very special configurations. 

The most important property of the BA function is that it is an eigenfunction of the 
multidimensional algebraically integrable Schrodinger operator L , which in our case has 
the form 

\ ^ mairria + l)(a, a) 
(a, x) 



(see ii). 

When 2t is a Coxeter configuration, i.e. 21 consists of the reflection hyperplanes for 
some finite reflection group W with -invariant multiplicities, then the corresponding 
operator L is the Hamiltonian of the generalised quantum Calogero-Moser problem (after 
Olshanetsky and Perelomov |]^) with special integer- valued parameters. The existence of 
the BA function in this case was proved in with the help of Heckman's result 0. 

At that time it was believed that the Coxeter case is the only one, when if) does exist, 
but it turned out to be not the case. The first non-Coxeter examples have been found by 
the authors in [H (see also 

According to the general procedure proposed by Berest and Veselov in [|I^ this led 
to the new examples of the hyperbolic equations satisfying the Huygens' Principle in 
Hadamard's sense. Motivated by these results Berest and Lutsenko started the investi- 
gation of the case when the potential depends on two coordinates only and found other 



new examples of the huygensian equations [11|. Later Berest proved [O] that they have 



actually found all such equations under assumption that the potential is homogeneous of 
degree (—2) . Since a generic Berest-Lutsenko potential could not be described by the 
construction [^], this was the reason for us to revise it. 

In the Section 1 we give such a revised definition of the BA function, which can be 
derived from the corresponding Schrodinger equation and therefore covers all the possible 
cases. It is remarkable that there exists an effective way to check for a given configuration 
whether BA function exists or not. Namely, as we prove in Sections 2 and 3, the following 
overdetermined system of algebraic equations is a necessary and sufficient condition for 
the existence of the Baker-Akhiezer function: 

mii(mf + l)(l3,B)(a,Bf'-^ 



E 



(/3,x)2i+i 



on the hyperplane (a,x) = (2) 
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for each a & A and j = 1,2, . . . , rria ■ 

They are equivalent to the vanishing of the first rria odd terms in the Laurent expan- 
sion of the corresponding potential 

ma(ma + l){a,a) , , 

at the hyperplane {a,x) = . Similar characterisation of the rational finite-gap potentials 
in one dimension has been first proposed in the famous paper |13| by Airault, McKean 



and Moser, who introduced the term "locus" in this situation. We will also use this 
terminology, calling the equations @ as well as its general affine version (see below) 



as locus equations. Duistermaat and Griinbaum discovered the interpretation of 
such equations as a trivial monodromy condition for the corresponding one-dimensional 
Schrodinger equation in the complex domain. We give a similar interpretation for our 
locus equations (0) in Section 2. 

However, to describe all the configurations, satisfying the locus equations (H) (the locus 
configurations) seems to be a very difficult problem. At the moment it is solved only in 
dimension 2, where the answer is given by the Berest-Lutsenko construction. In dimension 
n > 2 all known examples of the locus configurations are the Coxeter configurations and 
their special "deformations" In the section 4 we present all the results, which are 

known in this direction so far. 

The generalisation of our construction to the affine configurations of the hyperplanes 
is discussed in the section 5. The potential u and the locus equations in that case have 
the form: 

mi{mi + l){ai,ai) 
mj{mj + V){aj, aj){ai, aj)^^"^ _ ^ 



((a„a:) + c,)^^+i 



(5) 



identically on the hyperplane (a^, a;) + q = for all i = 1, . . . , and s = 1, . . . , . 
Unfortunately, so far a little is known about the affine locus configurations, which are 
not linear, i.e. with not all the hyperplanes to pass through one point. Apart from the 
one-dimensional case investigated in |13[, [|15], there are only some reducible examples 



discovered by Berest and Winternitz [|1^. In fact, we show that the classification problem 
for the affine locus configurations can be reduced to the linear case (0) by the isotropic 
projectivisation procedure. 

In the last section we discuss the relations of our BA function ip and locus configura- 
tions to the Huygens' Principle. The main result says that for any locus configuration in 
dimension n the corresponding hyperbolic equation 

(□at+i + m(xi, . . . = (6) 

satisfies Huygens' Principle for large enough odd N . Conversely, we show that if the 
equation (^ satisfies Huygens' Principle and all the Hadamard's coefficients are rational 
functions, then u{x) has a form (Hf) for some locus configuration. 
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We conjecture that this construction gives all huygensian equations of the form (□^+1+ 
u{xi, . . . , Xn))4> = . In the case n = 1 it is well-known result by Stellmacher and Lagnese 
[0. When n = 2 and u is homogeneous this follows from Berest's theorem The 
proof of the general case would lead to the solution of the famous Hadamard's problem 
in the class @. 

1. Rational Baker- Akhiezer function related to a con- 
figuration of hyperplanes. 

Let ^ be a finite set of noncoUinear vectors a = (ai, . . . , G C" with multiplicities 
TTia G N . We will assume that (a, a) = J2^=i (^1 ^ ■ 

Definition. A function ip{k,x), /c,x G will he called Baker- Akhiezer function 
(BA function), if the following two conditions are fulfilled: 

1) ip{k,x) has a form 

^(fc,^) = ^^|^e('=-), (7) 

where A{k) = Yla^A (^iO^)"*" ; P{k,x) is a polynomial in k with the highest term A{k) ; 

2) for all a e A 

d^{tlj{k,x){k,ar-) = dl{tlj{k,x){k,a)^-) = ... = dl'^--\xlj{k,x){k,ar-)=0 (8) 

on the hyperplane Ha '■ {k, a) = , where da = (a,-^) is the normal derivative for this 
hyperplane. 

Notice that (|^ means that ip is a rational function of k with the prescribed poles 
along the hyperplanes IIq,, a E A and with the asymptotic behaviour at infinity: 

^ = (l + o(l))e('='^) 

when /c — > 00 along the rays outside the singularities (cf.[|^). 

First of all, in the same way as in one can prove the following 

Theorem 1.1. If the Baker- Akhiezer function ip exists then it is unique and satisfies 

the algebraically integrable Schrddinger equation 

Lt^ = (9) 

where 

I- A I V ^"("^« + !)(«.«) 

Algebraic integrability of the operator (0) means that L is a part of a rich (supercom- 
plete) commutative ring of partial differential operators (see for precise definitions). 
This ring is described by the following theorem. 

Theorem 1.2. Let TZjs, he the ring of polynomials f{k) satisfying the following 
properties 

dj{k) = dif{k) = ... = d':^--'f{k) ^ (11) 
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on the hyperplane {a,k) = for any a & A . 

If the Baker-Akhiezer function ip{k,x) exists then for any polynomial f{k) G TZj^ 
there exists some differential operator Lf{x, such that 

Lfipik^x) = f{k)ip{k,x). 

All such operators form a commutative ring isomorphic to the ring TZ_a ■ The Schrddinger 
operator ( jJ^ j corresponds to f{k) = —k"^ . 

We give the proof of these statements in a more general affine situation in the Section 

5. 

We should note that there exists the following explicit formula for Lf (due to Yu.Be- 
rest ig). 

Theorem 1.3. The commuting partial differential operators Lf for f G TZj^ are 

given by the formula 



Lf = CN{adL)''[f{x)l (12) 

where = {—1)^/2^N\ , N = degf , f is the operator of multiplication by f{x) , and 
[adi)^ means the N -th iteration of the standard ad -procedure, adAB = AB — BA. 

The proof follows from the results of the next section (see Corollary 2.5). 

We should note that originally in another axiomatics for the - function was 
proposed. There was considered a function 0(fc, x) of the form 

0(fc,x) = P(A;,x)e('='^), (13) 

where P{k,x), as in (0), is a polynomial in k with the highest term A{k) , with the 
property 

x)) = dim, x)) = ... = dl^--\m x)) ^ (14) 

at the hyperplane 11^ . 

Comparing ([T3|), ([T^ with (0), (||) we see that the difference between these two ax- 
iomatics is due to the additional factor H/g^a (fc,/?)™'' . In the Coxeter situation consid- 
ered in [Q (see Section 4 below) this factor is not essential because of its symmetry. 

It turns out that this minor change makes the axiomatics less restrictive and leads to 
a richer class of the integrable Schrodinger operators. We will prove (see Corollary 2.7) 
that if there exists cf) satisfying the conditions ([TB|), (|T^) then there exists also the BA 



function ip with the properties (|^), (||) and in that case ip = -^j^ . Converse is not true: 
there are configurations, for which ip does exist but (p does not (see remark after the 
proof of the Theorem 4.4). 



2. Monodromy and BA functions. 

Let L = —A + u{x) be a Schrodinger operator with a meromorphic potential u{x) having 
a pole along the hyperplane IIc^ : (a, x) = , which is assumed to be non-isotropic: 
(a, a) 7^ . 



5 



We are looking for a formal solution of the Schrodinger equation L(j) = X(f) in the 
form 

(jy{x) = j2<i'^^\^^^r'"^ (15) 

s>0 

for some fi , where the coefficients ^i""* = (j)^f\x-^) are some analytic functions on the 
hyperplane , x-^ is orthogonal projection of x onto Ua , (j)'^^ ^ . 

Let's suppose that the equation Lc/) = has a solution of the form (|15D with some 
H < . Then the substitution into the equation gives immediately that the potential u{x) 
must have a second order pole along IIq : the Laurent expansion in the normal direction 
a has the form 

k>-2 

with cfl = — l){a, a) . 

Moreover, we obtain the following recurrent relations for the coefficients ^i""* : 

s-2 

(a, a)(/i(/i -1) - + s){fi + s - l))(f)s = (A + A)0,_2 - ^ Ci0,_i_2, (17) 

i=-l 

(s = 1,2,... ), where A is the Laplacian A restricted to the hyperplane IIq, and we 
omitted all the indices a in the coefficients. 

If 2/i ^ Z we can determine all (ps from (^) and obtain the solution (|TB|) starting 
from an arbitrary function (pQ (the same procedure gives also another solution with /i' = 
-1-/X). 

In the one-dimensional case this is a classical way (going back to Frobenius, see e.g. 
[p7[| ) to construct the basis of solutions of the corresponding equation 

—ip" + u{x)ip = Xip (18) 

in a vicinity of its regular singular point. In the case when the equation ( |I8D has no 
monodromy in the complex domain, i.e. all the solutions are single-valued, we have that 

1) fj, must be an integer: /i = — m, m G Z+ , 

2) the first 2m + 1 equations from ( |r7| ) must be compatible. 

In case if this is true for each energy level A we will say that the Schrodinger operator 
has trivial monodromy. 

In the multidimensional case there exists a generalisation of Frobenius's theory for the 
partial differential equations with the regular singularities in the complex domain (see 
[^). For the Schrodinger equation with a singularity along a hypersurface the regularity 
condition means that the potential has a second order pole at most. 

The considerations above motivate the following 

Definition. We say that a Schrodinger operator L = —A + u{x) with meromorphic 
potential u{x) having a second order pole along the hyperplane IIq, : {a,x) = has local 
trivial monodromy around this hyperplane if 
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1) the Laurent coefficient c^f^l in the expansion has the form c^f^\ = ma{ma + 
1)(q;, a) for some rria G Z+ , 

2) the system (jl^ with fi = —nia is compatible for any function 0o for all 

A e C. 

Theorem 2.1. L has local trivial monodromy around Ila if and only if the coefficients 
of the normal Laurent expansion of the potential u{x) near n„ 

s>-2 

satisfy the following conditions: c_2 = ma{ma + l)(a, a) for some G , and 

= = 4"^ = . . . = 41_i = on n„. (19) 

In that case the Laurent expansions of the corresponding eigenf unctions ^TB^ ) satisfy 
the conditions 

= 0?^ = • • • = 0?i-i = on n„. (20) 



Proof is similar to the one- dimensional case considered by J.Duistermaat and A.Griin- 
baum |T^. Let's demonstrate the idea in the simplest case when = 1 . After substi- 
tution (^) into the Schrodinger equation, we deduce that /x = 2 and derive the following 



recurrent relations for 



(-2 + C„2)0O = 

201 + c_i0o = 

202 + (-A - A)0o + co0o + c_i0i = 

003 + (-A - A)0i + Ci0o + Co01 + C_i( 



(21) 



where A is the Laplacian A restricted to the hyperplane 11 (we omitted all the subindices 
a in these formulas and assumed that (a, a) = 1 ). These relations allow one to find all 
the coefficients uniquely except 0o (which is an arbitrary function) and 03 , provided the 
consistency of the first four equations. From the first equation it follows that c_2 = 2 . 
Expressing 0i and 02 from the second and the third equations and substituting them 
into the fourth one we arrive at the relation 



~ 1 1 ~ 1 

;-A - A)(--c_i0o) - 2^-i(-A - A)0o + (ci - coc_i + -cl^) 



0, 



which should be valid for all 0o and A . Vanishing of the leading term in A gives 
c„i0o = , i.e. c_i = . The relation reduces after that to Ci0o = , thus ci = . 
Notice that as it follows from the second equation 0i = — ic_i02 = 0. This completes 
the proof in the case when ma = 1 . In the general case one should use induction 
e 111, p.196). 

One can consider a more general case, when u{x) has a singularity along an 
arbitrary hypersurface v?(x) = . However, analysis of the corresponding relations (pT|) 
shows that the hypersurface has to be a hyperplane (cf. ||I9|). 



arguments 
Remark. 
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Now let's consider a Schrodinger operator (|I]) L , corresponding to some Baker - 
Akhiezer function if) . We claim that such an operator has local trivial monodromy around 
all the singular hyperplanes. To prove one can consider for a given A the (n — 1) - 
dimensional family of the solutions of the Schrodinger equation 

(L - A)(/? = 

of the form (p = ilj{k,x) with k"^ = —A. They have proper pole behaviour near the 
hyperplane {a,x) = 0. Unfortunately, ijj^"^ depends on k and is not an arbitrary 
function on the hyperplane, so we have to present additional arguments. We'll prove a 
slightly more general result, which we will use also in the section 6. 

Theorem 2.2. Let the Schrodinger operator L = —A + u{x) have an eigenfunction 
ip{k, x) 

= -k'^ip 

of the form ip = P{k^x)e^^'^^ , where P is a finite sum of some functions which are ho- 
mogeneous in k and meromorphic in x . Then the singularities of u{x) are second order 
poles located on a union of non-isotropic hyperplanes and L has local trivial monodromy 
around these hyperplanes. 

Proof. The fact that singularities of u{x) must be located on the hyperplanes was 
proved by Yu.Yu. Berest and A. P. Veselov in ||20[ under assumption that P is a polyno- 
mial in k , but their proof works also in the case when P is a finite sum of the homogeneous 
in k functions. The fact that these hyperplanes must be non-isotropic follows from the 



zero-residue lemma of the same paper (see also [0). 

Let's now prove that the conditions (|T9|) are to be satisfied. After a proper choice 
of orthonormal basis we may assume that the hyperplane under consideration has the 
equation xi = , and let's consider the Laurent expansion for the function ip{k,x) : 

+ 00 

%l){k, x) = xr™ ^^^^^ a;2, . . . , Xn)x\. (22) 

1=0 

Let's prove first that m has to be positive. Let P" be the highest homogeneous term of 
P , then from the Schrodinger equation we have ^ kid/dxiP^ = . So P^{k, x + kt) is 
constant while t varies, hence if P" vanishes on the hyperplane Xi = , then it vanishes 
identically. Thus, P" and therefore ijj can not be zero at the hyperplane, so m in (p^) 
must be positive. 

Substitution (^) to Schrodinger equation immediately gives that c_2 = m(m + 1) 
and leads to the following recurrence relations: 

^ j 

(m(m + 1) - (j + 2 - m)(j + 1 - m))^,+2 = (A - k^)ijj - ^ di^^.i, (23) 

i=-l 

( j = -1, 0, 1, 2, ... ), A = + . . . -I- . To prove (|19D let's suppose that c_i = Ci = 



8x2'^ ' ' ' dXr, 

. . . = C2p_3 = , but C2p_i 7^ for some p < m + 1 . Considering j = —1, 1,3,... , 2p — 3 
it is easy to see that tpi = tp^ = . . . = tp2p-i = . From the form of the function 
Ip it follows that tjjj = Pj{k, X2, ■ ■ ■ ,Xn)e^'''^\ where Pj are finite sum of homoge- 
neous functions in k , k = (^2, . . . , kn) , x = (x2, . . . , x„) . Let P? be the highest 
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homogeneous term of Pj . By induction one can prove that P^j = {-lyk^^^Pl^aj and 

P2j-\ = {—^y~^kl^^~^~^^PoC2p-ibj , where the constant aj > and bi = b2 = . . . = bp = 
(by assumption) and bj > for m > j > p + 1 . Indeed, for Pjj follows easily from 
the relations (|23D. For -P2j_i o'^^ can use induction arguments similar to [1^ ( prop. 3.3, 



p. 196). 

Now let's consider the equation pi) with the resonance value j = 2m — 1 



(A - k^)tp2m-i - Ci^: 



2m- 1 

2m— 1— i- 



Since this holds identically for all k the highest homogeneous term should vanish. Simple 
calculation shows that this term is equal to 

~{P2m-l^l + -P2m-2pC2p-l) = ("l)" ^^^k^^"^ ^^PQ°(6m + am)c2p-l- 

Since 6^ + c^m > and Pq 7^ it vanishes only if C2p-i = . This completes the proof. 

It is remarkable that the BA function turns out to be symmetric with respect to k 
and X . For Coxeter configurations this property has been established in |^ . 

Theorem 2.3. Baker- Akhiezer function ip{k,x) is symmetric with respect to x and 
k : ip{k, x) = ip{x, k) . 

Proof. The idea is to show that tp{x,k) is also the BA function and then to use 
the uniqueness (theorem 1.1). Let's prove that ^^^^jl'^'^^ is a polynomial in x with the 
highest term A{x) , where A{x) and P{k,x) are the same as in (^. For that let us 
consider the conditions (|^) for ip{k, x) . They give a linear system for the coefficients 
of the polynomial P with the coefficients, which are proportional to the degrees (a, x) , 
a & A. Since this system has a unique solution, these coefficients are rational in x . 
Let' denote by Pj{k,x) the homogeneous term of ^-^j^ of degree —j in k. In terms of 
Pj{k, x) one can rewrite the equation (^ in the following recurrent way 

" d 

LP,{k,x) = 2j2h—P,+,, Po(fc,x) = l. 

From this it follows by induction that all the singularities of ip{k,x) in x belong to our 
configuration of the hyperplanes (a, x) = . Analyzing Laurent expansions for u{x) and 
Tpik^x) on these hyperplanes we conclude that ipik^x) has a pole of order rria along the 
hyperplanes {a,x) = 0. All that means that A{x)P{k,x) is a polynomial in x. But 
from the uniqueness of BA-function it follows easily that Pj{k, x) is also homogeneous in 
X with the same degree —j . Hence the highest term in x of the polynomial A{x)P{k,x) 
is equal to A{x)A{k) . Thus ^{x, k) = "^^(^ " g*-*^'"^-* . Properties of the Laurent expansions 
in X follow immediately from the theorems 2.1, 2.2. So we have all the conditions for 
ip{x, k) to be a BA function. The theorem is proved. 

Corollary 2.4. Baker- Akhiezer function ip satisfies the following bispectral problem 

Lix, -^Mk, x) = -e^k, x), Lik, x) = -x^ij{k, x), (24) 

where L is the Schrddinger operator ([7 
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Now we are able to prove Theorem 1.3. 

Corollary 2.5. The Baker-Akhiezer function ip is an eigenfunction of the operator 
for any / G . 

Proof. Due to the Theorem 1.2 and to the symmetry of ip for any / G TZ^x there 
exists a differential operator A(k, ^) such that A{k, ■§j:)ip = f{x)'ip . On the other hand, 
L(x, -§^)i^ = —k'^ip . Now we can use the identity (1.8) from [|I^ which states in that case 
that 

(adLYifM = i-adPyiAM 

for all r G Z+ . For r = N = oidA = degf the differential operator {—adk'^Y{A) in 
the right-hand side has zero order and is, in fact, the operator of multiplication by cf{k) 
with c = (— 2)^A^! . This means that ip is an eigenfunction of the operator {adLy'{f) 
with the eigenvalue cf{k) . This proves the theorem 1.3. 

Now let's explain why the existence of with the properties (|13|), (our old 
axiomatics, see section 1) implies the existence of BA function ip . This follows from the 
following general statement, showing that the new axiomatics is in some sense the most 
general one. 

Let A be any configuration of hyperplanes, L = —A + u{x) be a corresponding 
Schrodinger operator, A{k) = YlaeA^^^ k)™'" . Consider the functions (p of the form 

P is some polynomial in k and x: P = A{k)A{x) + . . . , where dots mean the terms of 
lower order both in k and in x . 

Theorem 2.6. If the Schrodinger equation Lip = —k'^ip has a solution (p of the form 
then (p{k,x) has to be BA function. 

Proof now is almost evident. Theorems 2.1 and 2.2 provide the conditions (^ for (p 
in a; -variable, and it has the required form (0) in x . Hence, (p{x, k) is BA function and 
according to the theorem 2.3 ip{x, k) = ip{k, x) . 

Corollary 2.7. If a function (p satisfies the conditions ^T^)-^L^) then ip = A~^{k)(p 
is the Baker-Akhiezer function l^-l^- 

Proof. As it follows from the results of the papers 0], [Q, the function (p must be 
an eigenfunction of the same equation (^. Then the arguments we used in the proof of 
the theorem 2.3 show that ip = A~^{k)(p satisfies the conditions of the theorem 2.6 and 
therefore is the Baker-Akhiezer function. 



3. Locus equations and the existence of BA function. 

Let A, as in Section 1, be a finite set of non-collinear vectors a G C" with given 
multiplicities ma G Z+ , 21 be the corresponding configurations of hyperplanes (a, k) = 
in C" and L = —A + u{x) be the Schrodinger operator with the potential 

u{x) = Y: "^"^7 + ^!^"'"^ (26) 

The theorems 2.1 and 2.2 from the previous section imply that if the BA function for 
the configuration 21 exists then in the normal Laurent expansions (|T6|) of the potential 
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u{x) the first odd terms c^'Li {j = 1, . . . yfria) should vanish identically on the hyper- 
plane {a,x) = . More explicitly, these conditions have the form of the following highly 
overdetermined algebraic system: 

Y: "^^"^^'^^^;f,f - on the hyperplane (a, x) = (27) 
for j = 1,2,... ,m„. 

We will call the equations (p7|) as locus equations, following Airault, McKean and 
Moser [|T^, who used this terminology in one- dimensional case. The configurations 2t 
which satisfy the locus equations we will call as locus configurations. 

The remarkable fact is that the locus equations (|27|) are not only necessary, but are 
also sufficient for the existence of the BA function. We will give the proof following the 
paper 



21 



Theorem 3.1. For any locus configuration 21 the BA function ijj{k,x) does exist 
and can be given by the following Berest's formula: 

^{k,x) = [{-2f^ M\A{k)]-\L + kY^[Y[{a,xr-exp{k,x)], (28) 

where M = ^^^^ m^, A{k) = Hae^l". ^)"'° ■ 

Proof. Let's consider the linear space V which consists of the functions 0(a;), x G 
C" , with the following analytic properties: 

1) ^i^)YlaeA(^^^)"^°' is holomorphic in C" ; 

2) for each a & A the Laurent expansion (|15]) for should not contain the terms of 
order —nia + 2 j — 1 ( j = 1, . . . , nia) , i.e. the conditions ( pO]) hold. 

The basic observation is the following 

Lemma. The space V defined above is invariant under the Schrddinger operator with 
the potential ^Mj ) provided that the locus conditions ^2li) are fulfilled. 



It follows easily from the imposed conditions on the Laurent expansions in a -direction 
for u{x) and (p E V . 

Now let's define the functions ipi {i = 0,1, .. .) in the following way: 



ifo = ]^(a;,a;)"'°exp(A;,x) 

a€A 



and 



y^i+i = {L + P)^.. (29) 

It's obvious that ipo belongs to V , hence by the lemma (pi also belongs to V . From the 
definition of these functions and the property 1 of it is clear that (pi can be presented 
in the form (pi = Ri{k,x)exp{k,x) , where Ri = Qi n«e^('^? ^)~™'" some polynomial 
Qi{k,x). From (p9|) it follows that the degrees of the polynomials Qi in x decrease: 
degQi+i < degQj . Therefore, for some (p^ ^ but v^at+i = (L + k'^)'pN = . Thus, 
(f) = ipN is an eigenfunction for the Schrodinger operator L. Let's prove that A^ in fact 
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equals to M = J2aeA • If we denote by the highest homogeneous terms of Ri in 
X , we see from (|^) that 

n 

From this we obtain immediately that for i = M = YlaeA 

= (-2)*^M! (30) 

From this we conclude that for i > M Ri (which is polynomial in k ) will be of the 
negative degree in x . Thus, it cannot be an eigenfunction for the Schrodinger operator 
L because of the following lemma due to F.A.Berezin p2[ . 

Lemma. If a quasipolynomial ip in k ip = P{k,x)exp{k,x) satisfies the Schrodinger 
equation (—A + u{x))il! = —k'^ip then the highest term in k of the polynomial P must 
be polynomial in x . 

This contradiction proves that the last non-zero function in the sequence (^) is 
(pM ■ Moreover, since (pM belongs to the space V we obtain using (|30| ) that ip{k,x) = 
{Rlj)~^ipM satisfies axiomatics (^,(H) in x as well as in k according to the theorem 2.3. 
So, we proved that ip{k,x) defined by the formula (^) is the BA function associated to 
a configuration A . 

Remark. The remarkable formula (^) for ip was discovered by Yu. Berest ([0), who 
proved that if %p does exist then it should have the form (p8|). 



4. Analysis of the locus equations and locus configu- 
rations. 

The next step would be to classify all the solutions of the locus equations {locus configu- 
rations). Unfortunately, this problem seems to be very difficult. In this section we present 
some results in this direction and all the known examples. 
4.1. Coxeter systems. 

The most natural examples of the locus configurations are given by the mirrors of the 
Coxeter groups. Recall that a Coxeter group W is by definition a finite group generated 
by some orthogonal refiections Sa{x) = x — with respect to hyperplanes in R" 

(see [^). If we consider all the refiections from the Coxeter group W , then the set 21 of 
the corresponding hyperplanes (a, x) = will be invariant under the action of W . The 
configuration 21 of these hyperplanes with arbitrary W -invariant multiplicities G 2+ 
gives an example of locus configuration. This fact follows immediately from the symmetry 
of the corresponding potential u{x) with respect to any refiection s^^cy. & A. 

In this case the Schrodinger operator L is the quantum Hamiltonian of the generalised 
Calogero - Moser system (see [0). The existence of the BA function for the root 



system of type An with ma = 1 was proved in [^], where some explicit formula for ip 
has been found. For the general Coxeter system it was done in using the Heckman's 
formula for the so-called shift operators in terms of the Dunkl operators Notice 
that our approach gives a new proof of this result. 



12 



Remark. In principle, one may try to extend these examples to the complex case, by 
considering a finite group generated by orthogonal refiections in complex Euclidean space. 
However, it is known (see e.g. ||2^) that all such groups are nothing but the complexified 



Coxeter groups. 

4.2. Deformed root systems. 

The first non-Coxeter locus configuration An{m) was introduced in [Q. It consists 
of the following vectors in R""*"^ : Cj — Cj with multiplicity m {I < i < j < n) and 
Cj — ^/men+l with multiplicity 1 {i = 1, . . . ,n). Notice that for m = 1 we have the 
root system An . We can allow the parameter m to be negative simply considering the 
vectors Cj — ej with the multiplicity —1 — m in that case (then, of course, we will have 
a complex configuration in C^~^^ ). 

Corresponding Schrodinger operator has the form: 



^ {Xi - XjY ^ {Xi - y/mXn+l)^ 

In the simplest nontrivial case n = 2 we have the following configuration (see figure 

!)• 




m 



, Q \ cosO = 

^ m + 1 



Figure 1 

The next example is related to the root system of C„ - type. Let's consider the 
following set of vectors in R""*"^ : 



C„+i(m,Z) 



Ci ± Cj with multiplicity k 

2ei with multiplicity m 

2^/ken+i with multiplicity I 

Ci ± y/kcn+i with multiplicity 1 
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where I and m are integer parameters such that k 



2m+l 
21+1 



e Z, I < i < j < n . In the 



case of 6*2(771, /) - system the parameters m, I can be arbitrary integers, the corresponding 
quantum problem was considered in ^. The corresponding configuration has the form 
shown on the figure 2. 



C2(m,l) 




Figure 2 

For n > 1 the corresponding Schrodinger operator has the form: 

m(m + 1) 



-A 



71+1 



^ 4k{k + l){x^ + x]) ^ 

^ {xf - X^jf ^ 



i<j 



i=l 



xt 



+ 



(32) 



1(1 + 1) ^j^4{k + l){x^ + kxl^,] 



n+l 



i=l 



(0^2 /^^2 \z 



where k 



2m+l 



2^^-^ . In the case / = m the system C„+i(m,/) coincides with the classical 
root system Cn+i (or -D„+i for / = m = 0). Again, as for the An{m) system, the 
parameters k,l,m may be negative, in that case the corresponding multiplicities in () 
should be —1 — k , — 1 — m or — 1 — / respectively. 

The simplest way to check the validity of the locus equations for these configurations 
is to use the following important property of the system (pTf): 

Theorem 4.1. A configuration 21 satisfies locus equations if and only if each 
two-dimensional subsystem of A gives a locus configuration. In other words, for each 
two-dimensional plane vr C C" the vectors a G ^ fl tt with their multiplicities nia must 
satisfy the locus equations. 

Remark. Notice the analogy with the similar property of the Coxeter and root systems. 

Proof. Let us denote by (3 the orthogonal projection of a vector (3 onto the hy- 
perplane (a, x) = , then (/3, x) = (/5, x) on this hyperplane. Let vr < a, 7 > denote 
the two-dimensional plane spanned by a and 7 7^ a . Then the subsum of (27) over 
/3 G TT < a,7 > becomes proportional to (7, x)"^-'"^ being restricted to the hyperplane 
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{a,x) = . All these subsums for different two-dimensional hyperplanes are independent, 
so we come to the following equivalent form of (pTf): 

for any two-dimensional plane vr G C" and for each a E ACm and j = 1, . . . , rria 

J2 mf^{m,3 + 1){(3, (3){a, (3y^-\(3,xy^^-^ = for {a, x) =0. (33) 

That gives the statement of the theorem. 

If we analyse the configurations An{m) and C„+i(m, /) from this point of view, we will 
have in each two-dimensional plane either a usual root system or one of their deformations 
A2{m) and C2{m,l) . For these two cases the locus equations can be checked by direct 
calculation. 

One can see that our configurations yl„(m) and Cn+i{m, I) have one common feature: 
they are obtained from Coxeter configurations by adding a special orbit of the Coxeter 
group with multiplicity 1 (a sort of "one-orbit deformation" of a Coxeter configuration). 
The following result demonstrates that such property is not accidental: the hyperplanes 
with large multiplicities always form a Coxeter subsystem. 

Definition. Let's say that the hyperplane 11/3 G 21 has a large multiplicity nip if in 
each two-dimensional plane containing vector f3 there are no more than mp + 1 vectors 
from A (without taking into account the multiplicities). 

Theorem 4.2. The set B C 2t o/ all hyperplanes with large multiplicities forms 
a Coxeter configuration and all other hyperplanes and their multiplicities are invariant 
under the action of this Coxeter group. 

Proof. We shall prove that for each j3 E B the corresponding reflection sp preserves 
the set 21 together with multiplicities. This implies, in particular, that sp{B) C B . To 
prove the invariance of 21 under sp let's consider as in theorem 4.1 an arbitrary two- 
dimensional plane TT, which contains /? , and the corresponding two-dimensional locus 
equation (p^): 

+ ^){l,l){l3,l?'-\l,x)-^^-\p,,)=^ = 0, 

where j = 1, . . . ,mp . Now we look at these equations for generic fixed linear 
system for unknowns = m^{m^ + 1)(7, 7)(/?, 7)(7, x)~^ of the form 

-(7^) l(A-)=o = 3 = 1,... ,mp. (34) 

We need the following elementary lemma: 

Lemma. // three unit vectors f3, 7, 7' belong to some two-dimensional subspace in 
and 

(/3,7)^ _ {f3,ir 
(7,x)2 (7',x)2 

for all X such that {(3,x) = then either 7 = ±7' or Sp{'j) = ±7' . 

Let's regroup the terms in (|34D into the groups corresponding to different values of 
l^'^ja . From the properties of the Vandermond determinant we easily conclude that the 
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sum of in each group should vanish. On the other hand, using lemma we see that 
there are only two terms in each group, and they correspond to the pairs of vectors 7, 7' 
with 5^3(7) = ±7'. Finally, we arrive to the condition + zy|(^ ,j.)=o = 0, which gives 
m^f{m^f + 1) = mY^my + 1) , i.e. = my . 

Remark. The Schrodinger operators (31), (^) remain integrable in a usual (Liouville) 
sense for the general (non-integer) values of the parameters m, / : there exists at least 
n = dimV^ independent commuting operators Li = L, L2, ... ,Ln- Indeed, for y4„(m) 
case (m is integer) it's easy to check that the polynomials ps = fc^ + A;| + . . . + fc^ + 
■m^k^_^_^{s = 1,2,...) satisfy the conditions ([TI|) and, according to the theorem 1.2, 
there exist differential operators Lg with the highest symbols Ps such that Ls^p = pgip 
and therefore [Ls,Lt] = (see the explicit formula (p!2D ). Since the coefficients of these 
operators depend on y/m in a rational way, one can define such operators for general m . 
For s = 2 one has the Schrodinger operator ([31|), and other Lg give its quantum integrals. 
In the case of C„+i(m, /)- system the similar arguments prove the integrability of the 
Schrodinger operator (p2D for the general l,m , and the commuting quantum integrals Lg 
have the symbols Ps = kj' + . . . + k^' + q'-^kl%^ {q = s = 1, 2, . . . ) . 

4.3. Locus configurations on the plane. 

Yu.Berest and I.Lutsenko [|l^ in the context of the Huygens' Principle have introduced 
the following family of the real potentials u on the real plane. In the polar coordinates 
they have a form 

2 d"^ 

u{r,<^) = - — -^logWlxM,. . . ,Xm{<^)], (35) 

where Xjiv^) = cos{kj{p + 9j), kM > . . . > > 0, A;^ G N, G R and • • • , Xm] is 

the Wronskian of Xi; • • • ; Xm ■ 

One can consider the natural complexification of the Berest-Lutsenko family in the 
following way. The set of all non-isotropic lines in is isomorphic to the cylinder 
C* ^ CP^\{0, 00} and can be parametrised by a complex parameter (^(modvr) 

X cos if + ysmip = 

Any configuration corresponds to a finite number of points in C* : ,93Ar with 

multiplicities mi, ... ,mN ■ The corresponding potential has the form 

--^E '^iP I (36) 

^ sm ((^ - yDj) 

where = + G C\{0} and (pijnodii) = arctan^. The complex Berest-Lutsenko 
potentials given by the formula (^Sj) with the complex parameters 9j , have the form (^) 
with ifj being the roots of the trigonometric polynomial W[(p] , their multiplicities are 
known to have a "triangular" form "^^^"J"*"^-' (see [0). 

Theorem 4.3. All the locus configurations on the plane are determined by the complex 
Berest-Lutsenko formula ([5^. 

Proof. First of all the locus equations (P7| ) in this case are equivalent to the following 
one-dimensional locus equations (cf. |T^) for the potential v{(p) = X]!Li ^n^T^^^^ ' 



d \ I X - mj{mj + 1) 
d^J V^sin2((^-(^^.) 



N, s = l,2,,... ,mi). 
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Now we can use the result from |2T|, which says that in its turn this is equivalent to 
the existence of the differential operator D with tt -periodic coefficients, intertwining the 
operator C = --^ + v{ip) with Cq = --^: 

CoD = DoCq. (37) 

The idea of the proof is close to the one demonstrated in the proof of the theorem 3.1, 
and we shall not reproduce it here. 

So, the only remaining thing to prove is that the relation ( p7|) implies that C can 
be obtained from Cq by classical Darboux transformations. Let's assume that D has 
the minimal order among all the intertwiners of C and Cq and consider its kernel: 
V = KerD . As it follows from ( p7D V is invariant under Cq : if Df = then D{Cof) = 
CDf = . Due to TT -periodicity of the coefficients of D , KerD is also invariant under 
the shift T : f{ip) /(<^ + 7r) . 

We would like to show that the spectrum of Co\v is simple and has the form 
(fcf , /c2, . . . , k\j) , where < fci < ^2 < . . . < fcjw are some integers. Suppose that there 
exists an eigenfunction f eV with the eigenvalue X k'^,k E Z. Since Cq commutes 
with T , we can assume that / is a Bloch eigenfunction: 

Cof = A/ 
Tf = fif 

If A 7^ fc^ / has to be pure exponent: / = Ce^^'^ or / = Ce"^^'^ . Since Df = 
the operator D can be factorised as 



D = DoF, 



dip f '' 



where D is a vr -periodic differential operator of the order one less than D (see e.g. 
When / = Ce^^^ we have F = ^± and CoDoF = DoFoCo = DoCqoF. 

Thus C o D = D o Co , so D is also an intertwiner with the order one less than the order 
of D. 

Thus the spectrum of £o|y consists only of the squares of integers: A = fc^, G Z . 
The same arguments show that A 7^ . So we have only to prove that the spectrum 
is simple. First of all there could be only one eigenfunction, corresponding to a given 
X = k"^ . Indeed, otherwise KerD contain the whole Ker{Co — A) and therefore D 
can be factorised D = Di o (£9 — A) with Di being another intertwiner of less order. 
Suppose that Cq has a Jordan block with X = k^ . Consider the Jordan basis /q, /i, . . . : 
(£0 ~ A)/o = 0, (£0 ~ A)/i = /o, . . . . Since /o can not be pure exponent (see above), 
/o = Acos{k(p + 9o) , then /i = ^ sm{k(p + 6o) + B cos{k(f + 6i) . Now from the invariance 
of KerD under the shift T we conclude that ^ sm{k(p + ^^o) also belongs to KerD . 
Together with /o the last function generates Ker{Co — A) , which leads to factorisation 
D = Di o [Co — A) and reducibility of D . 

Thus we have proven that KerD is generated by the functions xi,... ,Xn of the 
form Xj = cos{kjip + 6j) . The general formula (see e.g. |^) from the theory of Darboux 
transformations says that u = log • • • , Xn] ■ The theorem is proven. 
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We should mention that ahhough the formula ( P^ ) is explicit, it is not so easy to 
extract the geometric information about the locus configurations. For example, it is not 
clear how to prove the following theorem using this formula. 

It is very easy to show that all two-line locus configurations consist of two perpen- 
dicular lines with arbitrary multiplicities. Let's consider the first non-trivial case of 
three lines = 0, = and (7,0;) = 0, x G with arbitrary multiplicities 

ma,m^,m^ G Z_|_ , and ask when they form a locus configuration. Modulo the natural 
rotational equivalence we have the following classification. 

Theorem 4.4. All the three lines locus configurations are listed below: 

1 ) the Coxeter A2 configuration with multiplicities (m, m, m) ; 

2) the deformed A2{m) configuration with multiplicities (1, l,m) when m is positive 
and (1, 1, —m — 1) when m is negative; 

3) the three lines complex B crest- Lutsenko configurations, which can be parametrised in 
this case as: 

a = (1, a), /? = (1, 6), 7 = (0, 1) : - a6 + 6^ + 1 = 0, 

the multiplicities are (1, 1, 1) . 

Proof. Let 21 be an arbitrary three lines locus configuration. Let us consider the 
first case when 2t has at least two vectors with multiplicities greater than 1. Then the 
theorem 4.2 states that 21 has to be a Coxeter A2 - system. Now let us suppose that 
there is only one vector 7 = (0, 1) with multiplicity m > 1 . The theorem 4.2 states that 
other two vectors have to be symmetric with respect to the vector 7 , so we may fix the 



normalisation a = (1, A), (3 = (1, —A) . The locus equation (|27|) for a has a form: 



21 + A2 1-A2 m(m + 1)A 

^f^^ -+ ^7 =0 if x + Xy = 0. 

{x - Xyy 

From that it immediately follows that A can take only the following values: A = 
V2m+i ' V2m+i ' ^^^^ chcck that 2t is equivalent to the system A2{m) 

or A2(—m — 1) . The last case we have to consider is the case when all the three vectors 



a = (l,a),/5 = (1,6), 7 = (0, 1) have multiphcity 1. The locus equation (|27| ) for vector 7 
takes the form 

2a{a^ + 1) 26(62 + 1) 



{x + ay)^ {x + by) 



+ / J = if ?/ = 



or 



(a + 6)(a2 + 6^-06+1) = 0. 
The locus equations (pTl) corresponding to a and (3 can be written as follows: 

;i + a2)(l + a6) + 6(a- 6)3 = 
;i + 62)(1 + a6) + a(6-a)3 = 

In the case a + 6 = this system of equations is fullfilled if and only if = ^ , which 
implies that 21 is either Coxeter system A2 or deformed system A2{—2). In the case a'^ + 
b^ — ab-\-l = the above system holds automatically without any additional restrictions. 
Thus, the theorem is proven. 

Remark 1. We should mention that some of the configurations 3) contain an isotropic 
line (a = ±i,6 = or a = 0,6 = ±i) and therefore actually reduce to the two-line 
configurations. Notice also that when a = i/\/?> = —6 we have A2{—2) configuration. 
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Remark 2. It can be checked that for the configurations 3) from the theorem 4.4 the 
function with the properties (|T^ - |1^) doesn't exist. This demonstrates that the converse 
for the statement of the Corollary 2.7 is not true. 

Notice that from this result it follows that the locus of n lines is non-empty only for 
the special sets of multiplicities. Moreover, if the locus configuration is real than the set 
of multiplicities determines it uniquely up to rotation due to the following result. 

Theorem 4.5. There exists no more than one locus configuration in with given 
cyclically ordered set of multiplicities. 

Proof. Let A = {ai, . . . ,0^} be such a configuration for given set of multiplicities 
{nil, . . . ,mjv} , and let us fix normalisation = {— sin (pi, cos ipi), < (pi < < ■ ■ ■ < 
{p]\r < n . Considering the locus equations, we have, in particular, that 

sm'' [if j - <pi) 
Let's now introduce the function 



y-^ mi{mi + l)mj{mj + 1] 

^ Sm [Lfi - Lpj) 



We conclude that if $ = {ipi, ... , v^at) defines a locus configuration then necessarily 

d 



■U{ipi, . . . ,ipn) = 0. 



dip,, 

Function U being a sum of convex functions is a convex function in the domain < 
ipi < (P2 < ■ ■ ■ < 71 ■ Suppose it has one more extremum in the point $ = ((^i, ... ,(Pn)- 
Then U{ipi, . . . , ip^) should be a constant along the segment $ + ($ — < t < 1 , as 
well as each function (^j+i) _ prom that it follows that (pi = (pi + ipo for some 

Sin [ipi ^j) 

constant ip>o for all i . This means that system is defined uniquely up to a rotation. 

Corollary 4.6. // all the multiplicities are equal then the only real configuration on 
the plane is Coxeter, i.e. dihedral. 

The consideration of all two-dimensional subsystems implies the following more general 
result. 

Corollary 4.7 Any real locus configurations in R" with equal multiplicities must he 
Coxeter. 



5. Affine locus. 

In this section we present some results concerning the case, when the singular set of the 
potential u{x) of the Schrodinger operator is an affine configuration (3 of hyperplanes. 
So, we consider a Schrodinger operator L = —A + u{x) with rational potential having 
second order poles along some non-isotropic hyperplanes in C" . Let {as,x) + Cg = 
(s = 1, . . . ,K) be the equations of these hyperplanes. We will suppose also that the 
potential u{x) decays at infinity, i.e. u{x) while x — > oo along the rays outside 
singularities. 
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Impose now the condition that L has local trivial monodromy around its singularities. 
Then the Theorem 2.1 from the section 2 allows us to reformulate this condition as some 
algebraic conditions on the arrangement & of the singular hyperplanes {aj, x) + Cj = . 
First of all, it follows that the potential u{x) must be of the form 

^ mj{mj + l){aj,aj) 
u{x) = y ' — '-—^ (38) 

for some integers mi,... ,mx ■ Then the conditions (|T9|) imply that the Schrodinger 
operator with the potential of the form (|38| ) has local trivial monodromy around its 
singularities if and only if the following relations are satisfied: 



E 



((a„x) + c,)-+i ^''^ 

identically on the hyperplane (ctj, x) + q = for all i = 1, . . . ,K and s = 1, . . . , . 
We will call the relations ( pQ] ) as locus equations. The equations ( pTj) from the Section 



3 are their particular case, when all the hyperplanes pass through the origin. Sometimes 
we will refer to ( P^ ) and ( P7D as to affine and linear cases correspondently. 



As it follows from the section 2, the locus equations (p9D are necessary for the existence 
of a certain eigenfunction of the corresponding Schrodinger operator L (see theorem 2.2). 
As well as in the linear case (Section 3) the equations (|39| ) are sufficient for this. The 
following result has been proven in 

Theorem 5.1. Let L = ~A + u{x) be a Schrodinger operator with the potential of the 
form ( p^j which satisfies the affine locus equations ^3^). Then L has an eigenfunction 
(j) of the form (j){k,x) = P{k,x)exp{k,x) , where P is polynomial in k , Lcf) = —k'^ip. 

This eigenfunction (up to a normalization factor ) is given by the Berest's formula 
analoguous to (P^): 

K 

^{k,x) = [i-2)''M\C{k)]-\L + e'f[l[i{aj,x) + cjr^ exp{k,x)], (40) 

i=i 

where M = XljLi "^j ^ind C{k) = Y[f=ii'^ji ^)"^^ ■ The normalization is chosen in such 
way that ilj{k,x) = (1 + o{l))exp{k, clS 

We start the analysis of the affine locus equations and their solutions (locus configu- 
rations) from the one- dimensional case. 
5.1. One-dimensional case. 

In this case we have a configuration of K points zi, . . . ,zk with multiplicities 
mi, . . . , niK on the complex plane and the potential 



The locus equations in this case (for mj = 1 ) have been introduced in the paper by 
Airault, McKean and Moser [0. Duistermaat and Griinbaum [0] obtained them for 
the general multiplicities and proved that they are equivalent to the existence of the 
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differential operator D with rational coefficients, intertwining L = — 

T - ■ 

-^0 — —-j^ ■ 

Lo D = D o Lq. 



+ u(z) and 



All such operators L are the results of the classical Darboux transformations applied to 
Lq , so the potential u{z) can be given in this case in terms of the Wronskians by the 
well-known explicit formula: 

u{z) = -2— logiy[xi, . . . ,Xm], 

where the polynomials Xi? • • • ? Xm are defined by the recurrent relations Xi = ^7X2 — 
Xi, . . . , Xm = Xm-i (see Burchnall-Chaundy Adler-Moser |]l5l). The Wronskian is a 
polynomial Pm{z, ci, . . . , Cm) with the coefficients depending on the additional integration 
constants Ci, . . . ,Cm (see [|l5l for the details). 

Thus, the locus in the one-dimensional case is a union of the rational algebraic varieties 
of the dimensions m = 1,2,3,. 



parametrised by Ci, 



and the locus configu- 



rations are simply the roots of the corresponding Schur polynomials Pm{z,Ci, . . . , c^) • 
The solution ip of the corresponding Schrodinger equation —ip" + u{z)il) = —X'^ip has a 
form 



l + 5;]a,(z)A-Me"^ 



(41) 



i=l 



This is a degenerate rational case of the hyperelliptic BA function, corresponding to a 
general finite-gap operator 0. These rational BA functions (0) are characterized by 
the following properties in the spectral parameter (cf. [§3|). Let ^1, . . . be arbitrary 
parameters, ips be the Laurent coefficients of at A = : = X]^=^m ^^^si^) ■ Impose 
the following m linear conditions on the coefficients ip-m, ■ ■ ■ ? V'm-i '■ 



i'm-l + 6V'm-2s = 
s=l 
m—1 

m—2s—Z 

s=l 
m-2 

m— is— 4 

s=l 












(42) 



•m+l 







They are equivalent to a non-degenerate system for m unknown functions ai{z) and 
determine ip of the form (|l|) uniquely. The usual arguments 0], show that such 
a function satisfies the Schrodinger equation —ip" + u{z)ip = —X^ip with the rational 
potential 



u{z) = 2a[{z). 



(43) 



Notice that for given ^1,... the system (^) determines a m -dimensional linear 
subspace V{^i, . . . , C,m.) in C^"^ and therefore corresponds to a point of the Grassmannian 
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Gr(m, 2m) . It is more convinient to identify the system of conditions (^2]) with a point 
of some infinite-dimensional Grassmannian Grj^^ (see for the details). Namely, let's 
consider the linear space C[[A]] of formal series in A , and let be a subspace of C[[A]] 
with the following properties: 

1) A™C[A] G W G A~'"C[A] where C[A] is the space of polynomials and both inclusions 
have the same codimension m ; 

2) X'^WgW. 

We will suppose that the number m = m{W) in 1) cannot be reduced. The set of all 
such subspaces for m = 0, 1, 2, . . . we will denote as Gr^ following |33l . 



It is easy to see that the subspace of C[[A]] consisting of all Laurent series ifj = 
J2t^-m ^^'^s which satisfy the conditions (|^) represent nothing but a general point of 

(2) 

Gtq . In these notations the one-dimensional BA function corresponding to W is the 



unique element ipw of the form (^l]) which Laurent expansion at A = belongs to W 
for each z . We will denote by uw the corresponding potential (^3]). 

These considerations suggest the following extension of the axiomatics (0-|D of the 
multidimensional BA function. 

5.2. Equipped configurations and BA functions. 

Let A be again a finite set of non-collinear vectors in C" . We will prescribe to each 
vector a G A a subspace VT^") G Gr^^^ , and denote the corresponding integer ■m{W''"^) 
as rria . We will call the corresponding set of hyperplanes IIq, : (a, A;) = with the 
prescribed subspaces ly as equipped configuration 21 . 

Definition. For a given equipped configuration 21 the function tp{k,x) is called the 
Baker- Akhiezer function if it satisfies the following two conditions: 

1) ip has the form 

, . i!|£)e.-., ,44, 

where A{k) = Yla^A^'^^ /c)™" , P is a polynomial in k with the highest term A{k) ; 

2) for each a & A the Laurent expansion of ip in k in a -direction calculated at any 
point of the hyperplane Ha must belong to W^"'^ . 

Here by the Laurent expansion of a meromorphic function F{k) in a -direction at a 
point ko we mean the Laurent expansion of the function /(A) = F{kQ + Xa) at A = . 

If for each subspace W^°'^ the corresponding parameters ^ in (^) are zeros, our 
definition reduces to the definition of the BA function from the Section 1. Now we will 
prove the analogues of the theorems 1.1, 1.2 for a general equipped configuration. 

Theorem 5.2. If for a given equipped configuration 2t there exists BA function ip 
then it is unique and satisfies the Schrddinger equation 



-A + ^(a, a)ua{{a, x)) | ip 



-A;V, (45) 



where Ua{z) = Uyy(a){z) are the one- dimensional potentials, corresponding to the sub- 
spaces 

Theorem 5.3. Let TZ be the ring of polynomials f{k) with the following properties: 
for each a & A and any point k^ G 11^ the polynomial fa,ko{^) = f{ko+Xa) preserves 
the space .■ fa^.^^"^ C P^^") . 
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If the Baker- Akhiezer function tp{k,x) exists then for any polynomial f{k) G TZ there 



exists some differential operator Lf{x, such that 



Lfij{k,x) = f{k)ij{k,x). 

All such operators form a commutative ring isomorphic to the ring TZ . The Schrddinger 
operator corresponds to f{k) = —k"^ . 

The proofs of the theorems above follow in a standard way (cf. [^) from the following 
two lemmas. 

Lemma 1. // some function ip of the form (without the restrictions on the 

highest term of the polynomial P ) satisfies the conditions 2 from the definition of the BA 
function then the highest term in P must he divisible by A{k) = YlaeAi'^^ ■ 

Lemma 2. The BA function corresponding to an equipped configuration 21 has the 
following asymptotic behaviour at infinity: 



ip{k, x) 



exp{k,x) (l + y2a[''\{a,x))p^ + o{k- 



where a^^\z) are the first coefficients in the corresponding functions (^Ip ifja = 
and o{k~^) means the rational function of k which degree is less than -1. 

To prove the lemmas, let's expand ip in Laurent series in (a, k) on the hyperplane 
(a, A;) = . For convenience we may suppose that {a, a) = 1 and choose orthonormal 
basis in k such that {a,k) = ki , the other coordinates /c2, • • • ,kn we shall denote by k . 
Then up to unessential factor exp{k2X2 + . . . + fc„x„) -function (HI) takes the form: 



^{k,x) = e''^^^ klas{k,x), (46) 



and the Laurent coefficients are rational functions of k with possible singularities at 
zeros of homogeneous polynomial A{k) = k^'^A{k)\i:-^=Q . Since the sum '^s>-ma ^i^i 
is the Laurent expansion for ^^j^ , the degrees in k of its coefficients as decrease at 
s — oo (by definition, deg| = degp — degg ). Now we restrict our attention to the terms 

kfag with the maximal degree of in A; . From the remark above it follows that we have 
a finite number of such terms, and if we extract the highest homogeneous part in k in 
each term, we obtain the following finite expression 

4^%k^,x) = e^^'^ J2 Ka'sCk,x), (47) 

s>—ma 

where a° is the highest term in and all the have the same degree in k . It is clear 
now that constructed in that way ip'^ must obey the same restrictions (^21). This implies, 
in particular, that the sum (^) contains at least one term with s > . The outcome 
is that if we expand P{k,x) in the series in ki , P = J2j>oKPji^) 5 then extract 
from this sum the terms with the maximal degree in k , the result must contain at least 
one term with j > ■ Now let's present P as a sum of homogeneous in ki, . . . ,kn 
components P = Po + Pi + . . . , and suppose that the highest term Pq is not divisible by 
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k^" . In this case some other term Pi must contain A;™" , but its degree in k is clearly 
less than the degree of the term coming from Pq . This contradiction proves the lemma 
1. 

Moreover, in the extreme case when Pq has the form Pq = k'^'^Qo with Qo|a;i=o 7^ 
the reduced ip -function (^) up to a factor coincides with the one-dimensional BA function 
( PD tp{ki,Xi) . It's easy to see that this factor is simply (5o|fci=o • 

In particular, this implies that the second homogeneous term Pi in P{k, x) for the 
BA function ijj satisfies the following condition: 

where ai is the first coefficient in the corresponding one-dimensional BA function (|41|). 
We obtained this formula under assumption that (a, a) = 1 , in general it looks as follows: 

[(a,fc)^-™"Pi](^^,)^^ = (a,a)ai((a,x))[(a,fc)-"'^Po](„,,)=o- (48) 



Taking into account the restrictions ( ^8]) for all the hyperplanes {a,k) = , and we 
obtain that if Pq = A{k) = YI^ga^^^ ^)'"" ' 



aeA 



P, = A{k)J2at\i<y,x))^^y (49) 



which proves the lemma 2. 

Let's consider now for a given equipped configuration 21 the corresponding Schrodin- 



ger operator (|45|) . It is clear that the potential has the form (|38D . The corresponding 
affine configuration of the hyperplanes & we will call dual to the equipped configuration 
21. Suppose that the corresponding BA function does exist, then from the theorem 2.2 
we conclude that the Schrodinger operator (^) has local trivial monodromy and hence 
satisfies the locus equations (|39D . In other words, the dual configuration 6 must be a 
locus configuration. We believe that the converse is true, that is, each locus configuration 
appears in such way for appropriate BA function. The part 2 of the Theorem 5.6 below 
shows that each locus configuration is dual to some equipped configuration. So, the only 
problem is to check that for the function defined by the formula (1^) the properties 2 
from the definition of the BA function hold. Unfortunately, we couldn't find a proof for 
this. We can only remark that for all known affine locus configurations it is true. 
5.3. Geometry of afRne locus. 

First of all, it is easy to check that the following operations preserve the locus equations 
and therefore allow to produce the locus configurations: 

1) motions of the complex Euclidean space C" ; 

2) extensions of the configurations in C" to C"^ , m > n , induced by an orthogonal 
projection C™- C" ; 

3) union of two configurations which are orthogonal to each other. 

At the moment all known examples of the affine locus configurations can be con- 
structed using these operations from one-dimensional affine and multidimensional linear 
locus configurations. 

In particular, this is true for the configurations, corresponding to the operators intro- 
duced by Yu.Berest and P.Winternitz |16|. Analysis of these examples, however, reveals 



one more geometrical way to produce the locus configurations. 
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Let & be any afiine configuration of hyperplanes in C" . Let's imbed C" in C""^^ 
in the following way: x = {xi, . . . , x„) — > (xi, . . . , Xn, 1, 0) . For any hyperplane 11 in 
C" let's define tlie fiyperplane 11 in C""*"^ as a linear span of 11 C C" C C""^^ and the 
isotropic vector e = (0, . . . , 0, 1, . If (a, x) + c = is the equation of 11 in C" then 
the corresponding equation of 11 will be (a, x) + c(x„+i + iXn+2) = . 

Corresponding configuration & in C"^^ we will call as isotropic projectivisation of 

6. 

Theorem 5.4. T/ie isotropic projectivisation of an affine locus configuration <B in 
is a linear locus configuration G in C"+^ . 

Proof. We shall check the first of the locus equations for © , the others can be checked 
in the same way. So, we need to prove that on a hyperplane (a^, x) + Cs{xn+i + iXn+2) = 
the following identity holds: 

where oij denotes the normal vector of the hyperplane IIj C C"^^ . If Ilj C has the 
the normal vector aj = (aj, . . . , a") , then aj is the vector (aj, . . . , a", Cj, icj) . /^From 
that we immediately see that {aj,aj) = {aj,aj) and {as,Oij) = («s, «j) • Now since 
A = Xn+i +iXn+2 7^ almost everywhere on the hyperplane {ug, x) +Cs(x„+i +iXn+2) = 
we come to the identity 

rrij [rrij + 1 ) (ctj , ) (a^ , aj 



^ (K-,x) + c,A)3 = 



for (a^, x) + c^A = . But this identity after the rescaling x — > Ax takes the form 

EmAm^ -\- l)ia4,a.i)ia^,a^) 
\, ^\'' '\\ ' =0 for a„x +c. = 0, 
.^^ ((a„x) + c,)3 

which is exactly the first locus equation for the configuration G . 

Example. Let 6 be a direct sum of three-point one-dimensional configurations with 
the corresponding potential 



UiXi , . . . , Xj^ 



6xf — 12r,x, 



Then after the isotropic projectivisation we obtain the locus configuration with the po- 
tential of the form (cf. ||T6|): 

A 6x| - 12rj(x„+i + iXn+2Yxj 



tilXi, 



I Xn+2) — ^ 



In order to obtain a more general Berest-Winternitz's potential [16[ 



_ >, 6x^ — 12rj(x„4.i + ix„+2 + '^ji^^j 

(x^ + r,-(x„+i + zx„+2 + Cj)'^)^ 
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we should shift the pairwise-orthogonal triples of hyperplanes 

1 

Xj + {xn+i + iXn+2) =0 (j = 1, . . . , n) 

by Cj in Xn+i ■ 

Remark. The BA function in this example can be obtained easily using the following 
general remark. If ipi = Ri{k,x)exp{k,x) {i = 1,2) are given by the formula (^D|) for 
two orthogonal locus configurations &i and 62 then the function ip = RiR2exp{k,x) 
will correspond to the locus configuration © = ©i |J ©2 • This is clear from the structure 
of the formula (|40|). 

Thus, iterating such geometric procedures one can construct many new affine locus 
configurations. However, all of them are degenerate in the following sense. Let V{&) be 
the linear space of the normals to all the hyperplanes in & . We call G degenerate if the 
restriction of the complex Euclidean form on V{&) is degenerate. 

For a degenerate affine configuration one can define the following isotropic reduction 
procedure, which is inverse to the isotropic projectivisation. 

Let K be the kernel of the restriction of the Euclidean form onto V{&) . Consider 
the orthogonal complement V-^ of V in C" and choose a subspace L such that 

V + = K ®L. 

By an isotropic reduction of the degenerate configuration & we shall mean the configu- 
ration & n {a + L} , where {a + L} is a shift of L by a generic vector a G . 

Theorem 5.5. An isotropic reduction of a degenerate locus configuration is a non- 
degenerate locus configuration. 

The proof is similar to the case of isotropic projectivisation. 

These results may be interpreted in two ways. First, we can say that any affine locus 
configuration is a result of the isotropic reduction of some (degenerate) linear configura- 
tion. So, the classification problem for affine locus configurations reduces to the linear 
case. On the other hand, as we have shown, to classify all locus configurations it is 
sufficient to consider non-degenerate configurations only. Moreover, we can consider irre- 
ducible configurations only, i.e. exclude the unions of orthogonal subconfigurations. At 
the moment all the known non-degenerate irreducible locus configurations are linear or 
one-dimensional. It may well be the only possible examples. 

The following general result clarifies the geometrical structure of affine locus configu- 
rations. 

Theorem 5.6. Any affine locus configuration & has the following properties: 

(1) for each point xq G C" the subset ^ & of the hyperplanes passing through 
xq form a linear locus configuration; 

(2) for each hyperplane 11 G E the subset ©(H) C A of the hyperplanes parallel to 11 
forms an extended one- dimensional locus configuration. 

Conversely, any affine configuration with properties (1), (2) belongs to the locus. 
Proof. (1) Let's consider the locus equations for some hyperplane Ilj : (««, x) +Ci = 
passing through xq : 

> — — ^77 ' = ^ ^ n^, 50 
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s = 1, . . . ,mi. 

Now take x = Xq + y , then x G Ilj iff (a^, y) = and we have the following relation: 



Eiiij yiiij -|- 1 J ytXj, Uijj \u.i^ Uij j X ^ 

(rv. n,\2s+l 2^ 



for all y such that (ofi, ?/) = . Since the second sum is regular at y = , the first sum 
should vanish on the hyperplane {cei,y) = 0. Thus, we obtain linear locus equation for 
the configuration &xo ■ 

(2) To prove the second property, let's divide all the hyperplanes which are non-parallel 
to n into the subgroups in the following way: 11' and 11" belong to the same group if 
and only if their intersection is contained in 11 . Then in each group the sum of the 



corresponding terms in (|50D should vanish due to the property (1). The remaining terms 
is exactly the locus equation for the set of parallel planes ©(H) . 
The converse statement now is clear. 

We conclude this section by some negative results about locus configurations in R" . 

Theorem 5.7. For any locus configuration in the real plane there exists a point all 
the lines pass through. 

Proof. First we note that parallel lines cannot appear in locus configurations in . 
Indeed, the subset of parallel lines according to the previous theorem must give a real 
solution for the one-dimensional locus equations, which is impossible. 

Now let's fix some terminology: by vertices we will mean the intersection points for 
the lines from the configuration and by a ray - any ray from the configuration with the 
origin at some vertex (some rays may contain other vertices). Let's choose an orientation 
on the plane. This allows us to determine the oriented angle ip{li, h) between the ordered 
pair of raysZi,Z2, which varies from — tt to tt . We need the following property of the 
locus configurations in R^ : 

Lemma 1. For each ray li from the locus configuration in there exists another 
ray I2 with the same vertex and acute angle between li and I2 : 

Q<^{h,h)<\- 

Similarly, there exists a ray I3 with the same vertex such that — ^ < <^(/i,/3) < . 

Proof of the lemma follows from the linear locus equations ( P7| ) for the lines passing 
through a given vertex: it's clear that the sign of each term in it depends only on the sign 
of the cotangent of the oriented angle between a and /? . 

Lemma 2. Let li and I2 he chosen as in Lemma 1. Then if li contains another 
vertex of the configuration, the same is true for I2 . 

The proof follows from simple geometrical considerations. 

Let's consider now any vertex and all the rays of our configuration outgoing from this 
vertex. As it easily follows from the lemmas we have only two possibilities: 

1) there are no other vertices on these rays or 

2) there is at least one more vertex on each ray. 

Since we have a finite number of vertices, we obtain immediately that our configuration 
has only one vertex. Theorem is proven. 
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The same is probably true in R" but at the moment we can prove this only in the 
special case when all the multiplicities are equal. 

Theorem 5.8. Any affine locus configuration in R" with equal multiplicities is a 
linear Coxeter configuration. 

Proof. It's sufficient to prove that the configuration must be symmetric with respect 
to each its hyperplane. Since the parallel hyperplanes cannot appear in a real locus 
configuration, the statement follows from the theorem 5.6 and the corollary 4.7. 



6. Locus configurations and Huygens' Principle. 

Let us consider a linear hyperbolic equation 

Cip{x) = 0, £ = Djv+i + u{x), (51) 

where is the D'Alembert operator, = ^ • 

We say after J.Hadamard that it satisfies Huygens' Principle (HP) if its funda- 
mental solution is located on the characteristic conoid, i.e. this solution vanishes in the 
conoid's complement. 

Hadamard found some criterion for HP to be satisfied in terms of the so-called Ha- 
damard's coefficients Uu{x,^). They are uniquely determined by the following system of 
equations 

- + ''^'^ = -o^(f^-i) (52) 



j=0 



and the conditions that Uq{x,^) = 1 and U^{x,^) are regular at x = C, . These coefficients 
are symmetric with respect to x and ^ : Uy{x,^) = f/^(^, x) (for the details see the book 

Hadamard proved that the equation (|5TD satisfies Huygens' Principle if and only if 



is odd and U,\r = for z/ > ^ , where F = {(x,0 : (xo - 6)' - Etii^i - 6)' = 0} 
is the characteristic conoid. For the case when potential u (and, as a corollary, all the 
Hadamard's coefficients f/^ ) does not depend on at least one of the coordinates (say, Xq ), 
the Hadamard's criterion is equivalent to the condition U n-i = . 

2 

We consider the Hadamard's problem of description of all huygensian equations of the 
form: 

(□at+i + u(a;i, . . . ,XAr))v3 = (53) 

In fact, in our case for any locus configuration in C" the corresponding potential will 
depend only on the first n coordinates : u = u{xi, . . . , Xn), n < N . 

It turns out that huygensian equations of the form (|53D are closely related to the locus 
configurations. For the linear locus configuration in C" the corresponding potential 

nix) = Y: "^"("^" + ^)("'"). (54) 

is homogeneous of degree -2. 
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Theorem 6.1. For any real potential u{xi, . . . , x„) related to a linear locus config- 
uration the hyperbolic equation ( [J^ j satisfies HP if N is odd and N > 2 J2aeA^(^ ■ 
In that case the fundamental solution can be expressed via BA function. 

Conversely, if the hyperbolic equation ( |5^ with homogeneous potential u{x) : u{Xx) = 
X~^u{x) satisfies HP and all the Hadamard's coefficients are rational functions, then the 
potential u{x) must have the form ( |5^ for some linear locus configuration. 

Proof. The proof of the first statement repeats the arguments of the paper [|T0[] , 



where this result has been proven in the Coxeter case. It is based on the following relation 
between BA function and Hadamard's coefficients. If we have the Baker- Akhiezer function 
ijj of the form ([^, we can present it in the form 

x) = (f/o(e, x) + U^{^,x) + ... + Um{^, x))e(«'^), (55) 

where Uq = l,U^{x,C,) is homogeneous of degree —u in ^ , M = degA{k) = ^Qg_4^a- 
Since ip is symmetric in ^ and x (theorem 2.3), Ui, has the same degree in x. From 
the Schrodinger equation (H) for ip , Lip = —^^ip > L = —A + u{x) , we obtain: 



d 

—U, + L[f/,_i] = {u =1,... ,M + 1 with Um+i = 0). 



Since Ui; are homogeneous in x this implies the relations (p2D , so 11^, coincide with the 
Hadamard's coefficients. Now since Um+i = the Hadamard's criterion guarantees HP if 
N > 2M + 3 . Notice that it gives also the explicit formula for the Hadamard's coefficients 
and the fundamental solution for (^) (see for the details [p!0| 



Conversely, from the chain ( ^21) for the Hadamard's coefficients Up{x,^) for the ho- 
mogeneous potential u it follows that are also homogeneous in x (and, therefore, in 

f/,(Ax,0 = X-'U.^x^i) = U,{x, AO- 

This can be proven by the same calculation as in lemma 1 from |12[ , where the case n = 2 
was considered. Let's now consider the function ip defined by the formula (^). Then, 
from the Hadamard chain (^2|) and homogeneity of it follows in the same way as above 
that %p satisfies the Schrodinger equation 

(-A^ + ?z(x))^ = -e'V'- 

Notice that the potential u{x) must be rational since all the Hadamard's coefficients are 
supposed to be rational. This follows from the first equation of the Hadamard's chain 
(^2]). Now using the theorems 2.1 and 2.2 and the fact that u{x) is homogeneous of 
degree (-2) we conclude that u{x) has the form ( ^41) for some locus configuration. 
Remark. In the case when n = 2 i.e. u = u{xi,X2) , a stronger result (namely, without 
the assumption that the Hadamard's coefficients are rational) follows from the results by 
Yu.Berest and I.Lutsenko |TT|, |T^. 



Now let's consider an arbitrary (affine) locus configuration (5 such that the corre- 
sponding potential u{x) given by the formula (|38| ) is real for real x . This is equivalent 
to the condition 6 = S where © is a natural complex conjugation of a configuration 
& . The following result generalises the Theorem 6.1 for the general (affine) locus config- 
urations. 
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Theorem 6.2. For any affine locus configuration G C C" with (3 = S the cor- 
responding hyperbolic equation ( [5^ j satisfies Huy gens' Principle if N is odd and large 



enough: N > 2M + 3, M = ruj . 

Conversely, if the equation ( |3^ satisfies Huy gens' Principle and all the Hadamard's 
coefficients are rational functions, then the potential u{x) must be of the form ^SBj) for 
some affine locus configuration. 

Proof. The first part of this theorem can be derived from the theorem 5.1 and the 



resuhs by Yu.Berest (see also |1T9[)- We would like, however, to present here another, 
more illuminating proof. It is based on a different idea which will help us to prove the 
second part also. The idea is to reduce the affine case to the linear one using the isotropic 
projectivisation procedure. 

The main observation is encapsulated in the following lemma. Let U^{x,C,) = 
0, 1, . . . ) be some analytic functions of 2n variables x = (xi, . . . , ^ = ■ ■ ■ ,^n) 
which satisfy the equations (|5^) with some potential u{x) . Let's define now the new 
functions depending on x = (xi, . . . , x„, Xn+i,Xn+2) and f = (^i, ... , ^^+2) : 

U^X, = i^n+l + tXn+2r''i^n+l + ^a+2)"'^f/.( —. , 7 ^^7^) (56) 



and 



u{x) = {Xn+l + iXn+2) ^- ). (57) 

Lemma. The relations for U^{x,S,) and u{x) are equivalent to the similar 
relations in x,^ for U^{x,^) and u{x) defined by the formulas ([5^ and ^5^). 
The proof is straightforward. 

Now suppose that we have the real potential u{x) related to some affine locus con- 
figuration © = (5 C C" . Then the potential u{x) defined by ( |57D corresponds to some 
locus configuration & C C""*"^ which is exactly the result of the isotropic projectivisation 
defined in the previous section (see Theorem 5.4). Thus, according to the Theorem 3.1 
the corresponding Schrodinger operator L = — A„+2 + 'u(x) in C"'*'^ has the BA function 



iIj{^,x) which is given by the formula (28). Therefore, ip can be presented in the form 
analogous to (|55D, 

ipilx) = mix) + Uiilx) + ... + UM{li))e^^''\ (58) 

where Uq = 1 and the components Ui,{x,^) are homogeneous of degree —u in ^ and 
X , non-singular for x = ^ and satisfy the relations (^) in x,C, with the potential u{x) . 
Now let's consider their restriction for Xn+i + iXn+2 = (,n+i + i(,n+2 = 1 , 



We claim that the formula ( pQ]) determines the Hadamard's coefficients for the initial 
potential u{x) . 

First of all, let's notice that this formula really determines some functions of x,^ 
only. This can be derived directly from the formula (p8[). Indeed, it's easy to see from 



the inductive procedure (E^) that the pre-exponent in the BA function (|28D is a linear 
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combination of the "monomial" terms Yl^eA^'^^^y^''^'^^^^'^"' "^i^^ some integers PaiQa- 
Thus, a;„+2, ^n+i5 will enter in ip only as combinations x„+i + ?x„+2 and ^n+i + 

iin+2 ■ This means that the coefficients defined by (|59| ) indeed do not depend on 
x„+2, ^n+i, • As a coroUary of the homogeneity of U^, in x and ^ we may 
invert the formula (^) and obtain that Uy are related to U,^ by the formula (|56D. Now 
using the lemma we get the equations (|52|) for . It is clear then from ( p9|) that t/o = 1 
and f/i, are non-singular when x = ^ . The last remark is that the procedure ( p9D gives 
us the real-valued functions t/j, of x,^ e in the case when the initial potential u{x) 
is real, © = (5 . 

So, for any affine locus configuration we constructed the Hadamard's coefficients 
for the corresponding hyperbolic equation (|^), and Um+i = . Applying the Hadamard's 
criterion, we obtain the first part of the theorem. 

To prove the inverse statement, we suppose that the hyperbolic equation (|53|) is huy- 
gensian and has rational Hadamard's coefficients U^, with Um+i = . In that case we can 
define the homogeneous functions Uy{x, ^) by the formula (^61). According to the lemma, 
they obey the equations (^2]) with the homogeneous potential (p7|). Then in the same 
way as in the Theorem 6.1, we conclude that the function (|58D satisfies the Schrodinger 
equation Lijj = —^'^ip with L = —An+2 + u{x) . Now using the Theorem 2.2 in the same 
way as in the theorem 6.1 we deduce that the potential u{x) must correspond to some 
(linear) locus configuration & of non-isotropic hyperplanes in C""*"^ . But in that case 
the initial potential u{x) (see the formula (|57D) will correspond to the isotropic reduction 
& of & which should satisfy the locus equations due to the theorem 5.5. The theorem 
is proven. 

Remark. We have assumed that the potential u of the hyperbolic equation does not 
depend on xq , but essentially we have used only the fact that the sequence of the 
Hadamard's coefficients terminates at some step M . Actually all the results of this 
section can be generalised formally for any equation of the form ( pT| ) (even with the com- 
plex potential), which possesses the last property. In that case the singularities of the 
potential should satisfy the locus equations in C^'^ with the complex Euclidean structure 
defined by the metrics diag(— 1, 1, . . . , 1) . 

We conjecture that any hyperbolic equation (Dat+i + u{x))ip = with terminating 
sequence of the Hadamard's coefficients has a rational potential u{x) which corresponds 
to some locus configuration in C^'^ . We have proved this under the assumption that 
the Hadamard's coefficients are rational. The proof of this conjecture would lead to the 
solution of the famous Hadamard's problem in the class (|53D. Until now this problem is 
solved only when u depends on one of the coordinates (K.Stellmacher, J.Lagnese 117]) 



and when u is homogeneous and depends on two of the coordinates (Yu.Berest |T^). 



7. Some other relations and generalisations. 

7.1. The Baker-Akhiezer function tlj{k,x) related to an equipped configuration has the 
following remarkable property: it satisfies a system of differential equations not only in 
X but also in k -variables. Corresponding bispectral property of the one-dimensional BA 
function has been observed in the fundamental paper by Duistermaat and Griinbaum P 
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Let ilj{k,x) be a BA function related to some equipped configuration 21, © be the 
corresponding dual configuration of the poles of the potential u{x) given by (|38D. 

Let TZ be the ring of polynomials defined in the Theorem 5.3. Define also the dual 
ring S as the ring of all polynomials q{x) in x , satisfying the relations 

[q{x)]\(as,x)+cs=0 = 



dx J 

for all j = 1,2,... , rUs and for all the hyperplanes of the configuration & . 

Theorem 7.1. For any p{k) G TZ and q{x) G S there exist the differential operators 
Lp{x,d/dx) and Mq{k,d/dk) such that the BA function ip{k,x) satisfies the following 
bispectral problem: 

Lp{x,d/dx)i){k,x) =p{k)'ilj{k,x) , . 

Mq{k, d/dk)ilj{k, x) = q{x)^{k, x) ^ ^ 

The existence of the operator Lp{x,d/dx) is claimed in the Theorem 5.3. The ex- 
istence of Mq{k,d/dk) follows from the characterisation of ip by its analytic properties 
in X . Namely, one can show that the BA function ilj{k,x) is the unique function of the 
form 

^ B{x) 

where B{x) = Y['I=iii^s, x) +0^)'"'' and the dots denote the polynomial in a; of a smaller 
degree, such that the following conditions are fulfiled: 



d \ 



for each j = 1,2,... ,ms and s = 1,... ,A^. The fact that the BA function satisfies 
these conditions follows from the Schrodinger equation (^51) and the theorem 2.2. 

7.2. Similar approach can be developed for the trigonometric versions of our Schrodinger 
operators (Q). As well as in the rational case discussed in the present paper, the axiomatics 
of has to be amended in order to cover the most general case. We intend to discuss such 
axiomatics in a separate paper. The corresponding locus conditions have been described 
in |21[]. The bispectral property for the corresponding BA functions results in difference 
operators in the spectral parameter, which can be viewed as deformations of the rational 
Ruijsenaars and Macdonald operators (see [pS] ). 

7.3. The most of the results of this paper can be generalised to the case when the potential 
u{x) of the Schrodinger operator is a matrix- valued function. The locus equations for 



that case in dimension 1 have been described in |^6[. Multidimensional case is considered 
in ^. 
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